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ABSTRACT 

Two-dimensional and quasJ-3D Navler-Stokes solvers have 
been used to predict the static and dynamic airload characteristics of 
airfoils. The following three turbulence models were used: (1) Baldwln- 
Lomax algebraic model, (2)Johnson-King ODE model for maximum 
turbulent shear stress and (3) A two equation k-e model with law-of- 
the-Uall boundary conditions. It was found that In attached flow the 
three models gave good agreement with experimental data. In 
unsteady separated flows, these models gave only a fair correlation 
with experimental data. 

INTRODUCTION 


The flow field surrounding modern rotorcraft and propeller 
configurations is highly complex, and Is dominated by three- 
dimensional effects, transonic flow, flow separation and unsteadiness, 
and can be properly modeled only through the numerical solution of 
the 3-D unsteady Navler-Stokes equations. Since full 3-D simulations 
are costly, historically researchers have used simpler 3-D analyses 
such as the lifting line theory, which use a table look up of 2-D steady 
and unsteady airfoil characteristics. The airfoil tables needed may 
come from carefully performed experiments, or from 2-D computer 
codes. To be useful, the 2-D computer codes should provide: 

1. Reliable prediction of airfoil static load data and dynamic 
stall characteristics. 

2. A method for evaluation of the flow yaw effects on airload 
characteristics. 

3. A suitable turbulence model for properly modeling 
separated flows. 


In this study, 2-D and quasl-3D computer codes have been 
developed capable of predicting the static and dynamic load 
characteristics of straight and swept wings. Three turbulence models 
are currently operational. These are: (1) Baldwln-Lomax algebraic 
model, (2) Johnson-Wng ODE model and (3) Two-equation k-< model. 
This paper describes the performance of the above two computer 
codes for a variety of steady and unsteady flow conditions. The effects 
of turbulence model on the predicted flow properties are aiso 
evaluated. 
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MATHEMATICAL AND NUMERICAL FORMULATION 

In this work, the unsteady 2-D and quasi-3D, Reynolds 
averaged, compressible Navler-Stokes equations are solved In a body- 
fitted coordinate system using an alternating direction Implicit (ADI) 
procedure. The mathematical formulation has been described in detail 
In References 1 and 2, and only a brief outline of the formulation Is 
given here. 

Governing Equations 

The equations governing unsteady three-dimensionaJ flow are 
tie full Navier-Stokes equations, and may be written In a Cartesian 
c oordinate system as: 


qj + F x + Gy + H 2 » R x + Sy + T z (1) 

Here, q Is the unknown flow properties vector; F, G and H are 
the inviscid flux vectors; R, S and T are the viscous terms. The purpose 
or: the present work is to compute the unsteady viscous flow over 
arbitrary configurations undergoing arbitrary motion. To facflitate this, 
the following general curvilinear coordinate system is used: 

£ -«(x.y.z,t) 

ij - <i (x,y,z.t) 

f - r (x.y,z,t) 

T - t 

( 2 ) 

In such a coordinate system, equation (1) may be written in 
the following strong conservation form: 

q, + ^ + % + H f - Rj + s , + \ (3) 

The quantities q , F , G , H , R , S , T are related to their 
Cartesian counterparts through the metrics of transformation. 

In such a general coordinate system, the airfoD surface maps 
or to a p * Constant surface. The radial direction Is associated with the 
n coordinate, and the direction normal to the airfoil maps onto £ * 
constant lines or planes. For a detailed description of the flow and flux 
vectors in the cartesian and transformed coordinate systems, the 
reader Is referred to Ref. 1 . 

M lflitS Sweep Assumption 

In many applications involving flow over rotor blades, and 
pr openers, the flow In regions away from the root and the tip may be 
assumed to be Invariant along the spanwise direction. That Is, a large 
spanwise or radial component of flow may exist, but the derivatives of 


1 



the flow properties along the radial or spanwise direction may be 
assumed negligible. This assumption is sometimes called the "Infinite 
sweep" assumption and Is often used In the aircraft Industry to 
compute three-dimensional boundary layers, and to investigate their 
stability characteristics. Under the assumptions of infinite sweep, the 
above equations become 

+ Hj. = + T^. (4) 


The equation set (4) consists of 5 equations, corresponding to 
the conservation of mass, energy, and u,v,w momentum along the x,y 
and z directions respectively. In the special case where the blade 
sweep angle is zero, and the yaw angle of the flow relative to the blade 
is zero, the spanwise or radial component of velocity v Is identically 
zero, and the momentum equation along the y- direction may be 
neglected, resulting in 4 equations. 

Solution Procedure 

The above equations are parabolic in time, and may be 
advanced in time using a suitable stable, dissipative scheme. In the 
present work, a formulation similar to that described by Steger [Ref. 
2] was used. Standard second order accurate central differences 
were used to approximate the spatial derivatives, and to compute 
the metrics of transformation. The highly non-linear flux terms F and 
H, which are unknown at a given time level ‘n’ were linearized about 
their values at a previous time level ‘n’. The time derivative was 
approximated as a first order accurate, two-point, backward 
difference. This leads to a system of simultaneous equations for the 
flow vector q n + 1 . These equations were re-expressed as a penta- 
diagonal matrix system of simultaneous equations for the ‘delta’ 
quantity (q n + 1 -q n )- The penta-diagonal equation system was 
approximately factored into a product of tridiagonal matrices using 
the Beam-Warming approximate factorization scheme, as discussed 
in Ref. 1. 

Artificial Viscosity Model 

The use of standard differences to approximate the spatial 
derivatives can give rise to growth of high frequency errors In the 
numerical solution with time. To control this growth, a set of artificial 
dissipation terms were added to the discretized equations. These 
dissipation terms used a combination of second and fourth order 
differences of the flow properties, In a manner discussed by Jameson 
et at. [Ref. 3]. 

TURBULENCE MODELS 

As stated earlier, three turbulence models were considered In 
this work. These models are briefly descried here. 

Baldwin-Lomax Model 

This model is patterned after the well known Cebecl-Smith 
model, and has been extensively used by a number of researchers 
[Ref. 2,4]. It uses a two layer formulation to model the eddy viscosity. 
In the inner layer, the following expression Is used. 

p t q I 2 | (uy • vx) | (5) 

where I Is the "mixing length" measured as the distance of the point 
from the nearest solid surface, modified by the classical van Driest 
damping term, and the von Karman's constant . In the outer layer, the 
eddy viscosity is written as 

“ I “ F max Y max < 6 > 

where and represent the turbulent velocity and 
length scales In tne outer part of the boundary layer. The quantity 
F max Is computed as the maximum of the following function: 


F(y) - y | (Uy-v x ) | [i-exp{-yp^ w /2qp w )J 71 

where y Is the distance of the point from the nearest solid wall. 
Y Is the y- location where F(y) reaches a rrvixlmum. The quantities 
p ™?and /i are the density, shear stress and viscosity respectively. The 
subscript W represents conditions at the wall. 

The KJebanoff Intermlttency factor Is used to drive the eddy 
viscosity to zero far away from the boundary layer. In the wake regions 
downstream of the blade trailing edge, Baldwin-Lomax model is used 
with minor modifications. For a detailed discussion of this model, the 
reader Is referred to Ref. 4. 

Johnson-Klna One Equation Model 

The Baldwin-Lomax model Is an equilibrium model In the 
sense that is assumes that the eddy viscosity Instantaneously adjusts 
to the local flow characteristics. The Baldwin-Lomax model thus does 
not take Into account the upstream eddy viscosity or turbulent kinetic 
energy values. Johnson- King model attempts to rectify this situation, 
by solving an ordinary differential equation (ODE) for the maximum 
turbulent kinetic energy within the boundary layer at a given 
streanrrwise location. This ODE may be thought of as a simplified form 
of the Reynolds stress equation .This ODE is solved as an Initial value 
problem by marching along the flow direction, and automatically 
brings Into the account the upstream history of the flow. 

The eddy viscosity is assumed to be 
-exp^ti/^to)!® 

where 

^•OAp C^y k m (9) 

and 

M t0 * 0.0160/? o(x) u 0 5 7 (10) 

Here D is the van Driest damping function, y Is the distance 
from the wall, and k Is the maximum Reynolds shear stress at the 
current x- location. Furthermore, 7 is the Webanoff Intermlttency 
factor, ue Is the edge velocity, $ * is the displacement thickness, and 
a(x) is a modeling parameter. 

An ordinary differential equation Is used to model the 
streamwise development of turbulent kinetic energy, and is given by 

k m 1/2 - ' D m L m / ("> 

The quantities D m , , a 1 are empirically prescribed as done 

In Reference 5. The subscript m over some of the quantities Indicates 
that these quantities are computed at the y- location where the local 
shear stress Is maximum. The quantity (k^)^. Is computed from 
equations 8 , 9 and 10 with the value of a(x]assumed unity. The 
quantity a (x) is Iteratively determined, so that the value of the eddy 
viscosity computed from equations 8 , 9 and 10 and the maximum 
shear stress computed using equation (11) satisfy the following 
relationship: 

" ^tm I p I ( u y + v x) I m ^ 

For an efficient Iterative procedure for computing a (x), and the 
empirical relationships used In the Johnson-KIng model, the reader is 
referred to Ref. 5. 

Gorskfsk-c Model 

The third turbulence model considered In this work Is the well 
known k-< model, Implemented with a set of wall boundary conditions 
proposed by Gorskl [Ref. 6]. This model requires numerical solution of 
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two partial differential equations for the Instantaneous values of 
turbulent kinetic energy k, and the dissipation rate < at every point In 
the flow field. These equations may be formally written as 

(pk) t + (puk) x + (pvk) y * (t^k^ + (u k ky) y + S 1 

M t + fcue) x + (pve) y * (u 0 « x ) x + (f* 0 * y ) y + S 2 (13) 

Here /i k and are eddy viscosities controlling the diffusion of 
k and t ; S 1 and S 2 are source terms which describe the production 
and dissipation rates of k and «. As in the case of the original flow 
equations, these equations and their three-dimensional counterparts 
may be cast in a strong conservation form In a moving, body-fitted, 
curvilinear coordinate system. In the present work these equations 
were solved as a 2 x 2 system of partial differential equations using an 
ADI procedure similar to that used to solve the mean flow equations, 
after the mean flow has been updated at a given time level. 

In the vicinity of the solid wall, the values of k and « computed 
from the above equations were overwritten with values computed from 
the following assumed relationships for k and e : 

k-Cy 2 
€ - Constant 

These constants were evaluated using the values of k and t 
computed at nodes well within the logarithmic region of the boundary 
layer, In a manner documented In detail by Gorskl [6J. 

RESULTS AND DISCUSSIONS 
Steady Flow Studies Using the Baldwln-Lomax Model 

A series of steady viscous calculations were performed using 
the Baldwin-Lomax model. In Figs. 1-5, a number of flow results are 
shown In the form of surface pressures, loads, skin friction and velocity 
profiles for a NACA 0012 airfoil and a supercritical RAE 2822 airfoil at a 
variety of flow conditions and compared with experimental data. The 
following flow conditions are considered: 

Case 1 :Surface pressure distribution over NACA 0012 airfoil at free 
stream Mach number 1^ ■ 0.301, Angle of Attack a = 13.5 degrees, 
Reynolds number based on chord, Re * 3.9 million. 

Case 2: Static load characteristics for a NACA 0012 airfoil at - 
0.301, Re * 3.9 million. 

Case 3: Predicted lift curve slope dCj/cfa versus freestream Mach 
number for a NACA00 12 airfoil at a - 0 deg., Re - 9.0 million. 

Case 4; Surface pressure distribution over a RAE 2822 airfoil at 
corrected - 0.73, correctedo - 2.79 deg., Re - 6.5 Million. 

Case 5: Skin friction coefficient and velocity profiles on the upper airfoil 
surface for a RAE 2822 airfoil at - 0.73, a - 2.79 deg., Re - 6.5 
million. 

As seen In Fig. 1, good agreement was found between the 
computed surface pressures and experiment [7] at a high angle of 
attack flow prior to static stall. The static load characteristics of the 
NACA 0012 airfoil shown in Fig. 2 indicate good prediction of static lift 
and moment stall. However, the drag was overpredicted at high 
angles of attack. In these calculations, the transition effects on loads 
prediction were simulated by enforcing the transition point at 0.05% 
chord on both the upper and lower surfaces, In Fig. 3, a good 
agreement for the lift curve slope is observed between the present 
predictions and the experiment [8J. For the super-critical RAE 
2822 airfoil it is found that the predicted surface pressures shown on 
figure 4 are In good agreement with the experiment [9] over most of 
the airfoil. The predicted pressure expansion near the leading edge on 
the upper airfoil surface Is not as strong as observed In the 


experiment, presumably because the flow was assumed to be 
turbulent over the entire airfoil surface In the present calculations. The 
rntegrated lift and drag from the experiment are 0.7433 and 0.0127, 
respectively, which compare well with the values of 0.7432 and 0.0134 
from the present code. Fig. 5 shows that the skin friction in predicted 
well except near the trailing edge. The computed velocity profiles at 
two chordwise locations are In good agreement with the experiment. 

Additional viscous flow calculations computed using this 
i olver may be found In Ref. 1. 

l iynamic Stall of Oscillating Airfoils Computed Using Baldwln-Lomax 

liodfi) 

A number of dynamic stall calculations have been performed 
for several 2-D airfoils and a swept wing configuration with infinitely 
long span. The airfoils analyzed In the 2-D mode are a NACA 0012 
a irfoil and two modem helicopter aJrfofl sections (the Sikorsky SC-1095 
a rfol and the Hughes HH-02 Airfoil). The 2-D dynamic stall results are 
presented in Figs. 6 Through 13. The angle of attack variation during 
the dynamic stall Is given by a * 15 * 10 cos(*>t). The free stream 
K'ach number In these three cases was 0.28, the Reynolds number 
based on chord was 3.7 Million, and the reduced frequency based on 
snml-chord was 0.15. 

The calculations were carried out for two cycles of airfoil 
pitching motion starting with a steady state solution at 5 degree angle 
of attack, to remove influence of flow transients. Here results for the 
second cycle are shown. In figure 6, the dynamic stall calculations are 
presented for two grids: a 157 x 58 grid and a 253 x 58 grid. 

As can be seen In Fig. 6, predicted toads using the coarse grid 
agree with experiments [7] reasonably well during the upstroke. 
Although the maximum lift coefficient Is underpredicted and moment 
prediction was less accurate, the theory still dearly captures the 
moment and lift stall. The leading edge separation occurred when 
the angle of attack Is 17 degree during the upstroke. A large "primary 
vortex formed subsequently. As the airfoil Incidence Increased, this 
primary vortex passed along the airfoil surface and was shed into the 
wake at an angle of attack around 25 degrees. During the downstroke, 
however, the agreement between the numerical results and 
experimental data is not good. A secondary vortex from the trailing 
edge was shed during the downstroke. It appears that the strength of 
the secondary vortex is overpredicted (which takes place at a between 
25 to 23 degrees). As the Incidence of the airfoil continues to 
decrease, the agreement between predicted and measured loads 
Improves. 


The dynamic stall predictions for the Hughes HH-02 airfoil and 
tf™ Sikorsky SC-1095 airfoil are shown In Figs. 7 and 8. Good 
ag reement between the theory and experiment is also observed during 
th*j upstroke for these two cases. During the downstroke, the 
computed results are only in fair agreement with the experimental 
data. Overall level of agreement for these two airfoils Is similar to that 
for the NACA 0012 airfoil. 

In order to evaluate the effects of yaw angle on static and 
dynamic load characteristics, a series of calculations which take into 
account the effect of the blade sweep relative to the freestream, and 
thus the effect of radial flow on the load characteristics, have been 
done. In Fig. 9, the static lift versus angle of attack are shown for a 
NACA 0012 airfoil at 0.3 Mach number and Re *2. 7 million at 30 
degree yaw angle. Comparisons with the experimental data of Carta 
(l(') are also given. Good agreement between the two sets of data is 
observed. 

In Figs. 10 and 11, the dynamic stall load characteristics of a 
NACA 0012 airfoil at 0 and 30 yaw angle are compared. The flow 
conditions are: -0.3, Re -2.7 million. For comparison, Carta’s 

results are also given. Only a qualitative agreement between the two 
seti of data is observed. The yaw effects on dynamic stall hysteresis 
loops are predicted to be less profound than observed in experiment. 
Several factors may contribute to the discrepancies. Firstly, the 
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location of transition Is somewhat different between the two sets of 
data. In the computation, the flow Is assumed to be fully turbulent over 
the entire airfoil. In Carta s experiments transition depended on factors 
such as surface roughness, mean flow turbulence level, etc. and Is 
difficult to model. Secondly, three-dimensional effects of the swept 
wing with finite span considered In the experiment may be Important 
and the present quasi three-dimensional approximation may not 
suitable for this configuration. 

TURBULENCE MODEL STUDIES 

In order to assess the effects of turbulence models on the 
prediction of separated flows, a number of turbulent flow solutions 
have been computed using the three turbulence models and have 
been compared with each other and also with available experimental 
data. Several steady flow situations were considered. These are: 
turbulent boundary layer over a flat plate; attached subsonic flow over 
airfoil; and separated transonic flow over an airfoil. Following these 
studies, calculations were made for an airfoil experiencing dynamic 
stall, using higher order turbulence models. 

Attached Flows 

First, some results are presented for a turbulent boundary 
layer developing over a flat plate at zero angle of attack. The flow is 
nearly Incompressible and the Reynolds number Is five million [11]. 
The calculations were performed on a computational region which 
extended half plate length upstream and downstream of the flat plate, 
and one plate length above the plate surface; 151 equally spacing 
nodes were used along the streamwise and 51 nodes were placed at 
geometrically increasing distance in the normal direction. The first 
point off the surface located within a y+ of 2. The predicted surface 
skin friction and velocity profile at mid-plate are shown In Fig. 12 for 
the three models. These three models predict results that compare 
very well with experiment. Note that velocity characteristics of the 
sublayer and buffer regions were correctly captured by these modeis. 
This resolution Is crucial to the success of turbulence model In 
resolving the near-wall turbulent flow characteristics. 

Next results are presented for a turbulent attached flow past 
an airfoil. Fig. 13 shows computed and experimental pressure 
distributions for the NACA 0012 airfoil at an attached flow condition 
(K^ * 0.301 , a - 13.5 degrees and Re = 3.9 million). The computed 
surface pressures are In good agreement with experimental data [7]. 
it should be mentioned that the Johnson-King model which requires 
solution of an ODE starting from a user-prescribed point has some 
uncertainty as to where it should be activated. The activating location 
should be a finite distance downstream of stagnation point or 
transition location. Several locations were tested (x/c ranging from 
0.08 to 0.25 ) and it was found the solution in this case was not 
affected by the starting location. 

Transonic Flows with Separation 

Two transonic flow cases exhibiting mild and strong 
separation have been computed and compared with de tailed turbulent 
flow measurements, in Fig. 14, the surface pressures are shown for a 
NACA 0012 airfoil at experimental conditions [8]:^ = 0.899, Angle of 
attack * 2.66 deg., Re = 9 Million. It is seen that the Baldwin-Lomax 
and the k-« models predict similar pressure distributions, with a shock 
predicted stronger than the measurements. The predicted pressure 
distribution using the Johnson-KIng model is dependent on where the 
non-equilibrium formulation is activated. The location x/c ■ 0.15 gives 
best agreement with measured surface pressures. However, It Is not 
clear how a starting location may be chosen. 

The NACA 64 A0 10 airfoil at a shock Induced stail condition 
[12] ( « 0.8, a - 6.2 deg., Re * 2 Million) has also been 

considered for detailed turbulent flow comparisons. As can be seen In 
Fig. 15, the predicted pressure distributions are not in good agreement 
with experiment for any of the models. All of them predicted too 
strong and aft a shock and too little pressure recovery. It was found 
that the Johnson-King model (the activation location of the non- 
equilibrium calculation was chosen at x/c = 0.15) predicted a shock 
location forward of the other models, it should be mentioned that all 
these models showed some unsteadiness (5%) in their loads 


predictions (buffeting). Mean velocity and Reynolds shear stress 
profiles for this case are compared In Fig. 16. Except In the close-to 
the-wall and near-wake regions, reasonable agreement wWi the 
measured mean-velocity profiles Is found for all the three models. 
Because of the thick boundary layer predicted and the underpredfctlon 
of pressure recovery, the predicted mean velocity profiles do not 
match well with experimental data close to the wall and this effect also 
extends to the near-wake region(x/c - 1.02). Poor predictions are 
found for Reynolds shear stress proffles for alt of the models. The 
shear stress peaks are underpredicted and their locations shift more 
closer to the wall compared to experimental data. 

The behavior of the k-« model Is somewhat similar to the 
Baldwin-Lomax algebraic model in this case. However, a thinner 
reversed flow is observed using the k-c model. Additionally, the k-< 
model shows a better prediction of the Reynolds shear stress profiles 
close to the wall. 

From the two transonic flow cases just described, it appears 
that the k-e model does not hold any noticeable advantage over the 
simple Baldwin-Lomax model In predicting the shock-induced 
separated flows. It has been pointed out by Lakshminarayana [13J 
that the k-e model does not predict the separation point or the 
reattachment any more accurately. Therefore, inaccurate predictions 
of velocity and Reynolds shear stress profiles in the separated regions 
were not surprising. 

As for the Johnson-King model, it is seen that the solutions 
depend on the choice of activation location of the non-equilibrium 
calculation at least In separated flows. This shows that the optimum 
location differs from flow to flow. For unsteady, separated flow such 
as the dynamic stall problem. It would be difficult to determine the 
optimum location during the oscillating airfoil motion. In addition, 
calculations using the Johnson-King model tend to show small 
oscillations about mean loads for the steady cases and require more 
Iterations than other two models to ensure convergence. Thus, It Is 
concluded that the Johnson-King model Is not suitable for predicting 
the unsteady, highly separated flows. In dynamic stall calculations to 
be discussed next, only the Baldwin-Lomax algebraic model and k-« 
equation model were used and compared. 

Unsteady. Highly Separated Row (Dynamic Stall Case) 

The dynamic stall calculation for NACA 0012 airfoil previously 
reported was repeated for an comparison of the k-e and the Baldwin- 
Lomax modeis for this complex flow. Predicted aerodynamic loads are 
presented in Fig. 17 and compared with experiment [7]. The k-e 
model predicts higher lift during the upstroke. During the downstroke, 
predictions using the k-e model show trends similar to the Baldwin- 
Lomax model, except that a smaller second vortex shedding (around 
24 deg.) and much stronger third vortex-shedding (around 15 deg.) 
were detected. Results from both models only show a qualitative 
agreement with the experiment during the downstroke. 

COMPUTER TIME REQUIREMENTS 

The calculations presented In this work were performed on a 
CRAY X/MP supercomputer at the NASA Lewis Research Center. The 
computer time for a viscous solution using the BaJdwin-Lomax model 
Is 0.28 second per time step using a 157 x 58 grid. A converged 
viscous steady solution requires 2000 to 3000 time steps when a 
space-varying time step technique Is employed (roughly equal to 850 
seconds). For dynamic stall cases, the computer time required for a 
full cycle on a 157x58 and a 253x58 grid was 4100 seconds and 6600 
seconds, respectively. The swept wing configuration, requires 20% 
more computing time than the 2-D version, because an additional 
(spanwise momentum) equation needs to be solved. Computer time 
for the Johnson-King model is about 0.30 seconds per time step, and 
requires about 500 to 1000 more time steps than the Baldwin-Lomax 
model to achieve a steady state. For the k-< model, one time step 
requires 0.34 seconds and the same number of time steps as the 
algebraic Baldwin-Lomax model are needed to achieve a steady 
solution. 
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CONCLUSIONS 

An efficient solution procedure has been used to provide 
Improved prediction of complex flow phenomena associated with rotor 
flows. The two-dimensional, and quasi-three dimensional, 
compressible, full Navier-Stokes equations have been solved using an 
ADI scheme. Numerical results show that good prediction of static 
loads and dynamic stall hysteresis loops of rotor blade sections was 
feasible. Evaluation of three eddy viscosity models have been made. 
In attached flows the three turbulence models considered gave good 
correlation with experimental data. For strongly separated flows, eddy 
viscosity models available Including the k-c model are not adequate. 
No clear trend could be found favoring the use of higher order 
turbulence models In separated flows. 
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1 Surface Pressure Distribution over a NACA 0012 Airfoil at « 
0 301, a - 13.5 degrees, Reynolds Number » 3.9 Million. 
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2 Computed and Measured Static Loads over a NACA 00t2 Airfoil at 
K ^ m 0.3 , Re - 3.9 Million. 
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3. Predicted and Measured Lift Curve Slope vs. Mach Number for a 
NACA 0012 Airfoil at a * 0 degrees, Re = 9 Million. 



4. Surface Pressure Distribution over a RAE 2822 Airfoil at M * 0725, 
a = 2.92 degrees. Re = 6.5 Million. 
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Velocity Profiles and Skin Friction Distribution for Conditions shown 
i Figure 4. 
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9. Lift vs. Angle of Attack over a NACA 0012 AlrfoH at 30 degree yaw. 
Chordwise Mach Number 0.346, Re - 2.77 Million. 
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10. Computed and Measured Loads for a NACA 0012 Alrfoi at Deep 
Dynamic Stall Condition. ^ - 0.3. Re - 2.77 Million . Reduced 
Frequency 0.124 , Zero Yaw. 
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11. Computed and Measured Loads for a NACA 0012 Airfoil at Deep 
Dynamic Stall Condition. - 0.3, Re * 2.77 Million , Reduced 
Frequency 0.124 , 30 degree Yaw. 
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13. Row over a NACA 0012 Airfoil at 13.5 degree angle of attack; 0.3 
Mach Number; 3.9 Million Reynolds Number. 
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16. Velocity and Reynolds Stress Profiles for Conditions shown on 
Figure 15. 
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